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Введение. В статье рассматриваются вопросы расчета деревянных арок с учетом нелинейной зависимости меж- 
ду напряжениями и мгновенными деформациями, а также ползучести и геометрической нелинейности. В осно- 
ву положено интегральное уравнение вязкоупругопластической модели наследственного старения, изначально 
предложенное А. Г. Тамразяном [1] для описания нелинейной ползучести бетона. 

Материалы и методы. Мера ползучести принимается в соответствии с работой И. Е. Прокоповича и 
В. А. Зедгенидзе [2] в виде суммы экспоненциальных функций. Показан переход от интегральной формы закона 
ползучести к дифференциальной. Связь между напряжениями и мгновенными деформациями для древесины 
при сжатии определяется формулой Герстнера, при растяжении принимается упругая работа. Решение выпол- 
няется при помощи метода конечных элементов в сочетании с методом Ньютона-Рафсона и методом Эйлера по 
схеме ступенчатого увеличения нагрузки с корректировкой матрицы жесткости с учетом изменения координат 
узлов с последовательным вычислением дополнительных перемещений узлов, которые обусловлены невязкой 
сил. Предложенный подход для повышения точности определения деформаций ползучести на каждом шаге 
допускает вместо метода Эйлера использовать метод Рунге-Кутты четвертого порядка. 

Результаты исследования. На основе вариационного принципа Лагранжа получены выражения для матрицы 
жесткости и вектора дополнительных фиктивных нагрузок, обусловленных ползучестью. Разработанная авто- 
рами методика реализована в виде программы в среде МАТГАВ. Приведены примеры расчета для шарнирно 
опертых по концам параболических арок без промежуточного шарнира и с промежуточным шарниром в сере- 
дине пролета под действием равномерно распределенной нагрузки. Выполнено сравнение результатов, получа- 
емых в вязкоупругопластической и вязкоупругой постановке. Достоверность результатов подтверждена расче- 
том в упругой постановке в программном комплексе АМЗУЪ. 

Обсуждение и заключения. Для рассмотренных арок установлено, что даже при нагрузке, близкой к мгновен- 
ной критической, рост перемещений во времени носит ограниченный характер. Таким образом, характер их 
работы в условиях ползучести сильно отличается от характера деформирования сжатых стержней. 


Ключевые слова: ползучесть, деревянная арка, геометрическая нелинейность, вязкоупругопластичность, метод 
конечных элементов. 
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Введение. Древесина относится к материалам, проявляющим свои нелинейные свойства, как при 
кратковременных, так и длительных воздействиях. Большинство существующих реологических моделей дерева 
устанавливает его мгновенные свойства на основе закона Гука, однако для сжатой древесины характерна 
нелинейная диаграмма упругопластического типа [3]. Впервые вопросы совместного учета мгновенных 
упругопластических свойств древесины и ее ползучести стал исследовать К. П. Пятикрестовский. Указанные 
свойства он объединил на основе метода длительного модуля деформации [4—7]. В данном методе уравнение 
ползучести содержит время в явном виде, что существенно ограничивает его использование в случае сложных 
режимов нагружения и переменных нагрузок. В литературе имеется значительное количество публикаций по 
совместному учету ползучести и мгновенной нелинейности деформирования древесины при расчете сжатых 
стержней [8—14]. В частности, в работе [8] рассматривается вопрос применимости к описанию нелинейной 
ползучести древесины принципа суперпозиции Больцмана. В работе [9] для сжатых стержней проводится 
сопоставление теоретических расчетов с экспериментальными данными, а в статье [10] помимо сжимающей 
продольной нагрузки учитывается поперечная нагрузка. В работах [11-13] задача устойчивости решается для 
стержневых элементов в составе рамной конструкции. В статье [14] для сжатых деревянных элементов 
выводится длительная критическая сила с учетом нелинейной ползучести. 

Проблема расчета с учетом ползучести и мгновенной нелинейности деформирования актуальна не 
только для отдельно взятых стержней и рам, но и для таких стержневых систем, как арки. Целью настоящей 
работы является разработка методики расчета арочных конструкций с учетом нелинейных свойств материала 
при кратковременных и длительных воздействиях, а также геометрической нелинейности. 

Материалы и методы. В качестве соотношения, определяющего связь между напряжениями и 
деформациями, воспользуемся уравнением вязкоупругопластической модели наследственного старения: 


Л] °('| ЭС(ьт) 
в [Го () ат, (1 
Е (т дт 
Данное уравнение впервые было предложено в [1] для моделирования нелинейной ползучести бетона. 
Функция Хо) здесь устанавливает связь между напряжением и мгновенной деформацией, С(ёт) — мера 
ползучести. Диаграмма напряжения-деформации сжатой древесины при кратковременном нагружении хорошо 


аппроксимируется формулой Герстнера [15], имеющей вид: 
2 


Е 
С-В - дуб. (2) 
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Сжимающие напряжения в формулу (2) подставляются со знаком «+». При растяжении для дерева 
имеет место линейная диаграмма вплоть до разрушения. Выразив из (2) деформацию через напряжение, 


2 
в в | 
Е, г. 


получим: 


(3) 
На основе (3) функцию напряжений Ко) можно записать в виде: 
Г(с)= гн|1-|-5 
(4) 


Для меры ползучести древесины воспользуемся формулой, предложенной в работе В. А. Зедгенидзе и 
И. Е. Прокоповича [2]: 
(1) (6, +Ае")]1-Ве"9 |, (5) 
где Су= 2,87.10? МПа", Ау = 10,95.10? МПа", В, = 11 у=у, = 0,15 сут". 
Уравнение (1) можно представить в форме: 
(@] 


&=— 4+8, 6 
й (6) 


где Е = < __ секущий модуль, 5" = -[/[® (5) ОС (т) 


(с) ы [бд 


Для меры ползучести в виде (5) деформацию ползучести можно записать в виде суммы двух 





Ат — деформация ползучести. 


составляющих: 
=== +=, & = св, И @ (г) Е аты в, = Але (=) "а (7) 
Дифференцируя по времени у, получим выражения для ое роста каждой составляющей: 
а (89-5); © = АЛИ", ь 


Выведем на основе (6) связь между внутренними усилиями и деформациями для элемента при 








совместном действии продольной силы и изгибающего момента. Модуль упругости будем предполагать 
функцией от координаты у, которая по высоте сечения меняется от -Й/2 до Й/2. На основе гипотезы плоских 
сечений полную деформацию запишем в виде: 
Е==&+УХ, (9) 
“у 


где го — деформация среднего слоя, Хх = де — изменение кривизны. 


Подставим далее (9) в (6). Выразив с через &, получим: 
о=Е(у)(=-=)=Е(у)(&+ух-='). (10) 


Продольная сила и изгибающий момент в элементе связаны с напряжением через следующие 
интегральные зависимости: 


М = | суча: м= | ом. (11) 


Здесь А — площадь поперечного сечения стержня. 
Подставим (10) в (11) и преобразуем полученные равенства к матричной форме: 


[мм (12) 


р , | ь ЕА Е 
где М№ =|В (У)= А, М =|Е ( У) уё ЧА, [р] ар: =] — матрица приведенных жесткостей, которые 
А А 
определяются по формулам: 
ЕА= |Е(У) ) ЧА, Еб = 1Е(У) УаА, ЕГ = 1Е(У) ) УаА. (13) 


Решение задачи с учетом физической и геометрической нелинейности авторами будет выполнено при 
помощи метода конечных элементов. Используется стержневой элемент, приведенный на рис. |1. Осевая 
деформация при учете геометрической нелинейности представляет сумму линейной и нелинейной 
составляющей: 
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„ а 1 
& = & +5 = ри в (14) 





ыы 


Рис. 1. Стержневой конечный элемент 


Для получения системы уравнений МКЭ воспользуемся вариационным принципом Лагранжа. 
Потенциальная энергия деформации (ПЭД) записывается в виде: 


1 1 | 
П=5 | ве" = 5 ЕС (& +8 +ух-е ) аи. (15) 
| И 
где 5” — упругая деформация, равная разности между полной деформацией и деформацией ползучести. 
Запишем (15) как сумму четырех интегралов: 
п [265 (= +5 + УХ-= му + [Е в (в +=. + Ух—=°) 4У+ 
(16) 
+[ (5) У(& += + ух-=')аУ- [90 (ва + У—=°) и 
Первый интеграл в (16) запишется в виде: 
1 
[Е (вр + а + ухе ‘ай = [= Е(у)(& +ух-=°) Ах + 
0 А 
(17) 
+| 5 -5 У) 4Аах= 5 Мах- | 5-5”. КАах, 
[еее [емо 
где / — длина конечного элемента. 
Второй интеграл в (16) представим в следующей форме: 
1 
[Е (вр +а" +ух-=°) а = [= ‘ве [Е(У) ААах+ [0(& ) ау+ 
0 А | (1 8) 


+| 5 ° х|Е (у) уаяах- Га [Е У)=' аАак, 


Вторым слагаемым в правой части (18) можно пренебречь, учитывая его более высокий порядок 
малости. Третий интеграл в (16) записывается в виде: 


[Е У) УХ(& += +УХ—в ‘аи = [ЕО Иа, + кг) алк 
0 А 


(19) 
1 1 1 
+ [вх Ебах = [Мхак+ [вх Ебах. 
0 0 0 
Четвертый интеграл в (16) представим в о виде: 
[Е (=, +5" +УХ-= )ай = ЕС ) = ААах+ 
1 
+=, Е у)аАак- |(=). ву) [Е Е(›)алк= (20) 
0 | 0 А 


1 1 
= рома мк №4 -|(=) Е()аи. 
0 0 0 И 
Слагаемое | (ЕЕ ( у)аТ при минимизации по вектору узловых перемещений обратится в нуль. 
И 


Окончательно выражение для ПЭД примет вид: 


[ [ 1 [ [ [ 
— 5 [ Мах +2| ва’ ЕАак+ | Миах+2 | в хЕЗак- | хМ`ах-| в Мах+ (21) 
0 0 0 0 0 0 
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1 1 № м 1 
г 5" Мах) = 58 х} [[%- |. | + №, | а, 


где № =&ЕА+х Е5-М№М — средняя продольная сила в элементе, Хх, — среднее изменение кривизны 





элемента. 
С учетом (12) формула (21) примет вид: 


рт № 
П=- |= 02115:—25 „г |Аах+-М |ах, 22 
ВЫ [ое ем. 25 
где {= ХР. 

Для перемещений конечного элемента примем следующую аппроксимацию: 
—_И 





и(х) =и + вы: р (23) 
у(х) = +9, х+ 0х’ +0,х’; (24) 
ф(х) = | = а, +2а,х+3Зо,х”. (25) 

Коэффициенты полинома (24) определяются путем подстановки в (24) и (25) координат узлов: 
0 
000 У, 
(0 
оо 0 | (6) ы 
| 0, |= 
гг у, (26) 
0123 || 
(0 


Из (26) вектор {0} ={0 ©, ©, ©} выражается следующим образом: 


0 1 0 00 0 


Е 


—1 


[о 0 у 0 100 0 у 
О 1 0 Ф, 3 э 3 1 | Ф, 
Ос = . — НИЕ ИИ —__ —_— |* — Фф ( 5 27 
Пе о О И ы 
ОТ 21 3/7 ф, в 1 2 ||, 
т р р р ф, 


С учетом (27) формула (24) примет вид: 
ЕЙ. (28) 


где “= хх Хх]. 


Вектор {=} запишем в форме: 











— Щи обо 00 
(|, |=|_ 4 {9 Чи}=[в{и}. (29) 
а? 4 
После подстановки (29) и (28) в (22) ПЭД запишется в виде: 
гр, | № ь Теа" а 
п-5 у Га овации [8] а м. М. о (30) 


После минимизации функционала Лагранжа по вектору узловых перемещений приходим к системе 


уравнений следующего вида: 
([«]+[к, {0} =} +2, (31) 








ау а 
где [К] - [ву [2 ||В]4х — матрица жесткости, | К, зы №. ‚| Чх — геометрическая матрица 
„ м а 
[ м" 
жесткости, ЕП [В] вы й — вклад в вектор нагрузки деформаций ползучести, {Е} — вектор внешних 
0 


узловых сил. 
Решение физически и геометрически нелинейной задачи осуществляется при помощи метода Ньютона- 


Рафсона. Первым этапом производится расчет при #= 0. Приращение нагрузки выполняется квазистатически 
небольшими порциями с последовательным вычислением дополнительных перемещений узлов, которые 
обусловлены невязкой сил, корректировкой на каждом шаге касательного модуля упругости и координат узлов. 
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Далее временной интервал, на котором выполняется расчет, разбивается на конечное число шагов по 
времени. Расчет на ползучесть выполняется аналогично расчету на статическую нагрузку. Приращение 
деформаций ползучести на шаге #+/41 может быть определено с использованием метода Эйлера 








а И (32) 
и би 


Результаты исследования. Представленные уравнения и алгоритм расчета реализованы авторами в 
виде программы в среде МаЧаЬ. Был выполнен расчет шарнирно неподвижно опертой по концам арки 
параболического очертания под действием равномерно распределенной по длине нагрузки (рис.2) при 
следующих исходных данных: Ро = 1,48-10* МПа, А = 55 МПа, Г = 16м, } = 3,2 м. Поперечное сечение арки 
принималось прямоугольным с размерами В = 10 см, Й = 15 см. Арка разбивалась по длине на 40 конечных 
элементов, по высоте сечение делилось на 100 отрезков, количество шагов по времени принималось равным 
600, а по нагрузке — 200. Максимальное количество итераций на каждом шаге — 20. На рис. 3 показан график 
зависимости максимального прогиба от нагрузки. Штриховой линии соответствует расчет в упругой 
постановке. Для контроля правильности результатов решение упругой задачи также выполнялось в 
программном комплексе АМЗУЪЗ с использованием стержневых конечных элементов ВЕАМ 188. Количество 
конечных элементов принималось такое же, как и в Майа. Существенной разницы результатов выявлено не 
было. При расчете в упругой постановке резкий рост перемещений, соответствующий потере устойчивости, 
наблюдается при а = 10,5кН/м, а при учете мгновенной нелинейности деформирования — при 4 = 10кН/м. 


Ч 


ГУБЕ ЕЕ ЕЕ 
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Рис. 2. Расчетная схема конструкции: 


4 — распределённая нагрузка на арку, / — стрела подъёма арки, [, — пролет арки 
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Рис. 3. Зависимость максимального прогиба арки от нагрузки при кратковременном нагружении 
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Расчет с учетом ползучести для данной арки показал, что даже при нагрузках, достаточно близких к 
мгновенной критической, рост перемещений носит ограниченный характер. На рис. 4 приведен график роста во 
времени прогиба в середине пролета при нагрузке д = 8кН/м. Штриховой линии на данном графике 
соответствует решение в вязкоупругой постановке. Существенной разницы в результатах не выявлено. 

Также при указанных выше исходных данных был выполнен расчет трехшарнирной арки с 
промежуточным шарниром в середине пролета. В этом случае мгновенная критическая нагрузка оказалась 
существенно ниже. При расчете в упругой постановке она составила 4 кН/м, а с учетом мгновенной 
нелинейности деформирования — 3,3 кН/м. Как и в предыдущем примере, даже при нагрузке, достаточно 
близкой к мгновенной критической, ползучесть носит затухающий характер. Кривые изменения во времени 
максимального прогиба при а = ЗкН/м приведены на рис. 5. Обозначения такие же, как и на предыдущем 


графике. 
о 
22 
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Рис. 4. Изменение прогиба в середине пролета во времени при нагрузке 4 = 8кН/м 
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Рис. 5. Рост во времени максимального прогиба для арки с промежуточным шарниром в середине пролета при 4 =3 кН/м 
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Обсуждение и заключения. Полученные уравнения и разработанная методика являются 
универсальными и допускают применение произвольных зависимостей между напряжениями и мгновенными 
деформации, а также произвольных выражений для меры ползучести. Это позволяет производить расчеты 
конструкций, изготовленных не только из дерева, но и из любого другого материала. В результате 
проведенного анализа ползучести деревянных арок установлено, что в отличие от сжатых стержней, для них 
ползучесть носит ограниченный характер даже при нагрузках, близких к мгновенной критической. 
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